
NT@UW-11-19 



o 

D 

in 

o 

m 



Oh 
£^ 

(N 
> 

cn 
in 

od 
O 



Improving the Volume Dependence of Two-Body Binding 
Energies Calculated with Lattice QCD 

Zohreh Davoudi and Martin J. Savage^ 

^Department of Physics, University of Washington, Seattle, WA 98195-1560. 

(Dated: October 3, 2011 - 1:38) 

Abstract 

Volume modifications to the binding of two-body systems in large cubic volumes of extent L de- 
pend upon the total momentum and exponentially upon the ratio of L to the size of the boosted 
system. Recent work by Bour et al determined the momentum dependence of the leading volume 
modifications to nonrelativistic systems with periodic boundary conditions imposed on the single- 
particle wavefunctions, enabling them to numerically determine the scattering of such bound states 
using a low-energy effective field theory and Liischer's finite-volume method. The calculation of 
bound nuclear systems directly from QCD using Lattice QCD has begun, and it is important to 
reduce the systematic uncertainty introduced into such calculations by the finite spatial extent of 
the gauge- field configurations. We extend the work of Bour et al from nonrelativistic quantum 
mechanics to quantum field theory by generalizing the work of Liischer and of Gottlieb and Rum- 
mukainen to boosted two-body bound states. The volume modifications to binding energies can 
be exponentially reduced from O (^e~'^^ /Lj to O (^e~^'^^ / in nonrelativistic systems (where k is 
the binding momentum of the state) by forming particular combinations of the binding energies 
determined in the four lowest-lying boosted systems. Relativistic corrections to this combination, 
and others, that violate the exponential reduction are determined. An analysis of what can be 
expected from Lattice QCD calculations of the deuteron is performed, the results of which are 
representative of a generic loosely bound system. 
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I. INTRODUCTION 



One of the major goals of nuclear physics research is to determine the properties and in- 
teractions of nucleons and nuclei, and more generally hadrons, directly from the underlying 
theory of the strong interactions, quantum chromodynamics (QCD). The only known way 
to accomplish this is to numerically evaluate the QCD path-integral using Lattice QCD 
(LQCD) in which the space-time continuum is replaced by a finite-volume grid, and the 
integrals over the fields at each point in space-time are performed using Monte-Carlo. The 
systematic errors introduced by a finite space-time volume, Lr" x T, with a finite lattice spac- 
ing, 6, can be systematically removed by performing calculations in multiple volumes, with 
multiple lattice spacings and using the theoretical knowledge of the associated functional de- 
pendences. The computational resources that are required for LQCD calculations with light 
quark masses (m^) near their physical values, with small enough lattice spacings (6A^ <^ 1 
where is the scale of chiral symmetry breaking), and in large enough lattice volumes 
(m7rL> 27r) that permit reliable extrapolations and interpolations of phenomenologically 
important observables are now becoming available with the deployment of multi-peta-flop 
machines. 

While calculating (post-dieting) the masses of the lowest-lying hadrons to high precision 
is computationally demanding (even without strong isospin breaking and electromagnetism) , 
calculating the binding energies of systems composed of the lowest-lying hadrons (nuclei) 
to the same level of precision is significantly more demanding. For instance, calculating the 
mass of the proton at the percent level requires precision ~ ±10 MeV, while calculating 
the deuteron binding energy at the percent level requires precision of ~ ±20 keV. So 
while the same gauge-field configurations can be used for both calculations, the statistical 
precision required in the evaluation of the relevant correlation functions differs by orders of 
magnitude. Further, while the error in the calculation of a hadron mass introduced by the 
finite lattice volume scales as ~ e"'"'"^ for large volumes [1], the error in the calculation of 
a binding energy is largely dictated by the size of the bound state [H [21 13]. For instance, 
for a two-body bound state, the finite volume error scales as ~ e~'^°^ /L, where is the 



binding momentum of the state, B = — yyml — Kg ± yrn^ — — mi — , where B is 

the infinite-volume binding energy. For a loosely bound system such as the deuteron, the 
deviation between the ground-state energy calculated in a given LQCD calculation and the 
actual binding energy will be dominated by the size of the deuteron ^ and will, in general, 
be significant except in very large volumes. These volume modifications arise from the 
exclusion of single-particle momentum modes in the bound state wavefunction due to the 
periodic boundary conditions (BC) that are imposed on the quark and gluon fields in the 
spatial directions. 

During the last year, LQCD calculations have observed bound systems of baryons. Ev- 
idence for a bound H-dibaryon (a state with the quantum numbers of AA) was found in 
Uf = 2 ± 1 LQCD calculations with ~ 390 MeV [U E], and subsequent evidence was 
found in = 3 calculations with ~ 840 MeV [5]. Also, evidence was reported for ^He, 
^He [7| and the deuteron in quenched calculations with ~ 800 MeV. The infinite 
volume binding energy for each of these nuclei was determined by calculating the ground 
state energy in a number ensembles of gauge-fields with different volumes and then extrap- 
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dating to infinite volume. While the generation of quenched gauge-fields is inexpensive 
computationally compared to the generation of QCD gauge-fields, the results of quenched 
calculations cannot be used to reliably predict quantities in nature. 

An important observation that was recently made by Bour et al [H] in nonrelativistic 
systems is that the volume modifications depend upon the momentum of the bound state 
in the lattice volume, as moving bound states have different momentum modes excluded 
from their two-body wavefunction. The implication of this, when extended to quantum field 
theory, is that the unextrapolated binding energies of composite systems, and in particular 
light nuclei, calculated with LQCD will depend upon the total momentum of the system. 
Bour et al |9] were interested in calculating the scattering of bound states in a numerical 
evaluation of a low-energy effective field theory path integral, and found that the momentum 
dependent contribution to the two-body ground state energy had to be removed prior to 
using Liischer's method (TUl E] to determine the phase-shift. 

In this work, we extend the quantum field theory formalism established by Liischer [101 [H] 
and generalized to moving systems by Rummukainen and Gottlieb |12J, to determine the 
volume modifications of binding energies of bound systems composed of two spinless particles 
with s-wave interactions moving in a finite cubic volume, extending the results obtained by 
Bour et al [9] in nonrelativistic quantum mechanics. By forming combinations of the ground 
state energies of two-body systems with different lattice momenta, volume modifications can 
be exponentially suppressed in the nonrelativistic limit, and we determine the violations of 
this exponential suppression. In the case of the deuteron, the lowest energy-eigenvalue in the 
np system with lattice momenta of P = with |dp = 0, 1, 2, 3 allow for the infinite-volume 
deuteron binding energy to be determined with a volume modification of AB^'""''^^ 20 keV for 
L> 12 fm, orders of magnitude smaller than in the ground state energy of the |d| = system 
alone. From a practical standpoint, forming such linear combinations of boosted ground state 
energies requires significantly smaller computational resources than computing the ground 
state energies in multiple lattice volumes. The former can be accomplished with one set 
of quark propagators on one ensemble of gauge fields, while the later requires generating 
multiple ensembles of gauge fields and quark propagators on each. In some sense, these 
linear combinations of binding energies represent an exponential "volume-improvement" of 
the binding energy calculation on a given ensemble of gauge-fields. 

An analysis of possible determinations of the deuteron binding energy from a single lattice 
volume was carried out in Ref. [13]. For a large enough lattice volume, the energies of the 
lowest two levels with |d| = fall below the t-channel cut, and the effective range expansion 
(ERE) of pcot6 can be used. Truncating the ERE at the first two terms, which is known 
to be a good approximation for nucleon-nucleon scattering, allows the deuteron binding 
energy to be obtained by an interpolation. Of course, it is desirable to not make such an 
approximation in extracting the deuteron binding energy, but this would require more states 
below the t-channel cut and hence even larger lattice volumes. 

II. SPINLESS PARTICLES WITH S-WAVE INTERACTIONS 

For spinless particles of mass mi and m2 interacting in an s-wave, the scattering amplitude 
below the inelastic threshold is uniquely specified by the phase-shift 6o, and is proportional 
to l/(p cot 5o — ip) where p is the magnitude of the momentum of each particle in the center- 
of-momentum (CoM) frame. When these two particles are confined to a cubic volume of 
spatial extent L subject to periodic BC's, and with total momentum P = ^d, the energy- 
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eigenvalues of the system are the solutions to ^ 



q* cot5(g*) 



Z^\l-q'\Kmi,) , 



(1) 
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where n is a triplet of integers. A " * " denotes a quantity determined in the CoM frame, 
and for a system with P = ^d and total energy the energy in the CoM frame is 
s = E*"^ = i?^ — |Pp which defines the 7-factor (which depends explicitly on d), 7 = E/E*. 
The magnitude of the three- momentum in the CoM frame, <?*, is determined by E*^ = 
{\l q*'^ + ml + \J g*2 + m'^)'^ and the factor a that appears in eq. (2) is 
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and a tilde over any variable denotes scaling 
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In eq. 0, A 

by a factor of L/(27r), e.g. g* = g*L/(27r). In the case of equal masses (mi = m2) a 
and the expressions in eq. ([T]) and eq. ^ reduce to the known result for boosted systems 
of equal mass [T6l [T7] . Further developments are required in order to recover the results 
obtained for nonrelativistic systems by Bour et al from eq. ([T]), eq. ^ and eq. (|3]), as we 
now describe. 

Assuming the scattering amplitude admits a single bound state in infinite volume, the 
location of the lowest energy-eigenvalue in a finite cubic volume is dictated by the behavior 



of Z^^Ij for q*'^ < 0. It is clear from the form of Z^^]^.j in eq. (2) that there are no poles 
along the negative axis, and the Poisson resummation formula can be used to determine its 
asymptotic behavior at large — g*^. It is straightforward to show that for g*^ < 
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where 



7m = 7m|| -|- mj_ = (7 — 1 



m ■ d 



d + m 



(5) 



and the m are triplets of integers. By setting q* = iK, the eigenvalue equation in eq. ([T| 
becomes 



pcot 



17m 



(6) 



^ Exponentially suppressed corrections to this relation for tttt and NN scattering of the form e 
been determined in Ref. ^14^ and Ref. ^15J, respectively. 
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In the infinite volume limit F^'^\kL) = 0, and the eigenvalue equation becomes 

Pcot(5(p)|p^,,^ + Ko = , (7) 

which correctly reproduces the location of the pole in the S-matrix. While it is not required 
for this analysis, below the t-channel cut p cot 6{p) can be expanded in powers of the CoM 
energy pcotS{p) = — ^ + |rp^ + defining the ERE, where a is the scattering length and 
r is the effective range. For the lowest few d vectors, the finite volume functions F^^\kL) 
are (where we keep in mind that 7 depends upon |d|) 

^(0,0,0)/^^) = 6e-«^ + 6v^e-^'^^ + ^ e"^'^^ + 3 e^^.L _^ 
^(0,0,1) (,^) = 2(^2e-- + S^^e---) 



+ 2V2 e"^'^^ + 2 cos (27r«) 



e 



8 cos (27ra) /^^^ / ..l , cos (47ra) ^^^^ \ 



= 2( e-"^ + 2 cos (27ra) J — ^e"V ^'^^ 

V 7 + l 



7 V3 + 72 



kL 



+ ^ ( e"'^''-^ + cos (Ana) 
73 V 



27^ + 1 



+ I + 2cos(47ra) J^^e-V"^'^^ 1+ ... 



7^ + 1 



F^^'^'^VkL) = 6cos(27ra) e 

V 7 + 2 



+ 3V2 e-^'^^ + cos {Aira) 



3 -^|{272+l)«L 



272 + 1 



e 



2 /cos(27ra) /^^^^r , , 9 -^hP+»^L 



+ 3cos(4vr«)^^e-V^«^ + ... , (8) 

where the ellipses denotes terms that scale as ~ ^-V^i^qL jj^ higher. In the large volume 
limit where the F*^'^-' will give rise to K^^^ that are close to kq, the k'^'^^ can be determined 
in a perturbative solution to eq. ([6]). Introducing the dimensionless parameter A, writing 
= Ko + X^^l^ + X^nf^ + ... along with the substitution F(^)(/€(^)L) ^ \F^'^\k^'^^ L) and 
equating orders in A, leads to 

4"^ = f F^''\^oL) 



= ^ , (9) 

where is the residue of the bound-state pole in the S-matrix, and the higher k^^^ can be 
determined in a similar way. For a given boost- vector, k^^^ scales as k^i^^ ~ e~'^°^/L and ^2*^' 
scales as ^2*^^ ~ e~'^'^°^/L. Consequently, the contributions to k^*^^ that scale as ~ e~'^°^/L, 
~ g-v^«o-Ly'^ g^j^j ^ ^-V^KaL jj^ originate in and are of the forms given in eq. (8). 

It is clear from the explicit expressions for F'^'^\kL) given in eq. (pj) that linear combi- 
nations that provide universal cancellations of finite volume effects in binding energies do 
not exist in general. This is due to the appearances of both 7-factors and a-factors that 
explicitly depend upon d. However, the nonrelativistic (NR) limit where 7 = 7^^^) = 1 
and neglecting the binding energy compared to the rest mass of the constituent hadrons, 
a = a^^^^ = nil/ [mi + "^2)? allows for relations to be constructed for a given value of 
^(NR)_ Corrections to the relations can then be constructed as an expansion in 7 — 7^'^^) 
and a — a'^^^^ The case of equal masses, mi = m2, is special because a = | for any binding 
energy and not just in the limit where the binding energy is small compared to the rest 
masses. In the NR limit it is useful to write 

F'-'H.L) ^ Y. /i"""""'' ^ . (10) 

where the coefficients /j'*-'*-"' for the case of equal masses, rrii = m2 with a^^^^ = |, are 
given in Table [l] and the case where m2 = 2mi with a^^^^ = | is given in Table |ll| The 



TABLE I: The coefficients f^'"'^^ in eq. nlOh that determine the leading four finite-volume correc- 



tions to the two-body binding energy in the nonrelativistic limit for a system with mi = m2. 
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ratios of coefficients in column in Table II reproduce the quantity r(k, |) that are 

tabulated in Table 1 in Ref. [9j. 

It is also worth pointing out that the ground state energy of the unboosted, |d| = 0, 
equal-mass, a = |, system has leading and sub leading volume corrections that are three 
times larger than those of the |d| = 1 and |d| = \/2 systems. So while it does not constitute 
an exponential reduction in the volume modifications, the binding energy of the |d| = 1 and 
|d| = a/2 systems will be significantly closer to the infinite volume binding energy than that 
of the |d| = system. 
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TABLE II: The coefficients /j^^^^a^ in eq. lilOh that determine the leading four finite- volume cor- 
rections to the two-body binding energy in the nonrelativistic limit for a system with m2 = 2mi. 
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A. Volume-Improvement for Equal Mass Systems : a = \ 



2 



The equal mass systems are special, as mentioned previously, because a = | is independent 
of the binding energy of the system, a feature that is not present for mi ^ m2. Using the 
coefficients in Table [T] it is straightforward to construct relations that eliminate the leading 
and subleading orders of the volume modifications. We consider five relations: 
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1 ( 3^,(0,0,1) _ ^(0,0,0) ^) 



= ^0 + (l + «:o^)e-^°^ + OUV'OL^,^-^] , (11) 

where 7^ = 1 + r^^ldp with rj = From the forms of the volume expansions given in 

eq. (11), we see that eliminates the largest three volume contributions up to relativistic 
corrections, while and 7t*" eliminate the first two contributions, and and eliminate 
only the first. 
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III. NUMERICAL EXPLORATION OF THE DEUTERON BINDING ENERGY 



The deuteron is the simplest nucleus, comprised of a neutron and a proton. Its binding 
energy is B = 2.224644(34) MeV which corresponds to a binding momentum of kq ~ 
45.70 MeV (using the isospin averaged nucleon mass of M^r = 938.92 MeV). As it is a 
spin-1 system composed of two spin-| nucleons, its wavefunction is an admixture of s-wave 
and d-wave. The system is predominantly s-wave with a small admixture of d-wave induced 
by the tensor {L = S = 2) interaction. For the purposes of this analysis we will neglect 
the small d-wave admixture in the deuteron wavefunction and assume that the deuteron is 
entirely s-wave. The low-energy s-wave scattering of a neutron and proton with J'^ = 1+ is 
well described by effective range theory, where the ERE of p cot 6 converges rapidly with just 
the first two terms. The scattering length is known to be a^^^^^ = 5.425(1) fm, the effective 



range is r 



1.749(8) fm, and the shape parameter is anomalously small and neglected. 
The central values of these two parameters in the s-wave amplitude give rise to a deuteron 
binding energy of -B ~ 2.212 MeV with a corresponding hq of ~ 45.58 MeV, which are 
within ~ 0.5% of the actual deuteron binding parameters. 
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FIG. 1: The Z^o\l; q*^,0) functions (blue curves) , normalized by , for the deuteron in a cubic 

volume with L = 10 fm. The dashed (red) curves correspond to —\/—q*'^, asymptotic form of the 
function as q*'^ — )• — oo. The hyperfine splitting of the lowest poles in the |dp = 2,3 functions, 
determined by 7 — 1, is evident. 



The functions 



2 ryjd), 

7L0F^OO 



^1; g*^, 0) for |dp = 0, 1, 2, 3 are shown in fig. jlj along with their 

asymptotic form — ^— g*^ at large — g*^. It is clear that, of the four functions shown, the 
|d| = function will give rise to the largest volume modifications to the deuteron binding 
energy, as reflected in the deviation from its asymptotic form. This is because this function 
has a pole at g*^ = 0, while the |dp = 1, 2, 3 functions do not. It is also clear from fig. [l| 
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and magnified in fig. |2| that the |d| = 1 function exhibits the smallest deviations from 
its asymptotic form and the volume modifications to the deuteron binding are expected 
to be the smallest of the four considered. An interesting point to note from fig. |2] is that 
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FIG. 2: The ;^y^zJo^(l; q*'\ 0) functions for the deuteron with |dp = 0, 1, 2, 3 in a cubic volume 

-K^j^j- The dotted (red) curve corresponds to the asymptotic 



with L = 10 fm evaluated at q*"^ 
form of — K for L — )• oo. 



despite ^oo'*^'^^ having volume modifications that start at O {^e '^^/Lj, there are significant 
cancellations between all of the exponential contributions, leaving the function very close to 
its asymptotic value over most of the range of k. Figure |3] shows the ground state energy in 
the deuteron channel (negative of the binding energy) as a function of the spatial extent of the 
volume through numerical solution of eq. ([6]). Also shown in this figure are the contributions 
from the O (e^'^"^/L) volume modifications, and from the volume modifications up to and 



including O (c'^'^o^/l). 

Forming the linear combinations of the k^^^ given in eq. (11), the k*, from the exact 
numerical solutions to eq. ([6]), gives rise to the improved estimates of the deuteron binding 
energy that are shown in fig. |4| Surprisingly, there is little difference between the volume 
modifications improved to O (^~'^'^°^ / Lj and those improved to O (^e~^'^°^ /Lj for volumes 

with 10 fm< L< 20 fm. For L> 12 fm the K\ except for , provide estimates of the 
deuteron binding energy that are significantly closer to its actual binding energy than the 
ground state of any given d spectrum. The combination is closer to the infinite volume 



binding energy than one would expect. While it is improved to O [e 



'\f2KQL 



/L ) it appears to 



be better than any of the others that are improved to higher orders. However, this is true 
only at these "intermediate" volumes, while in the very large volumes the predicted hierarchy 
is, in fact, found. In this combination there is a subtle cancellation between different volume 
dependences in the range of volumes that are shown in fig. |4j 
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FIG. 3: The ground state energy in the deuteron channel. The blue, purple, brown and gray solid 
curves are the exact energy for each system with total momentum defined by |d| determined from 
eq. ([g]). The solid red lines is the infinite volume ground state energy. The dotted curves result 
from using the analytic forms for the volume modifications at O given in eq. p-lv, while 

the dashed curves are from the analytic forms up to and including O 
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FIG. 4: The exponentially-improved estimates of the ground state energy in the deuteron channel 



from the relations given in eq. (11). The gray dashed lines denote the ground state energies in 



the systems with a given d- vector that are shown as the solid curves in fig. [3| 
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IV. VOLUME-IMPROVED FITTING 



While it is important to form the exponentially volume-improved combinations of binding 
momenta, it may not be the method that is actually implemented in the analysis of the 
results of LQCD calculations. The existence of the relations shows that the volume mod- 
ifications in a prediction of the deuteron binding energy from one ensemble of gauge-field 
configurations can be exponentially reduced (in the NR-limit) with minimal additional com- 
putational resources. However, this reduction can also be accomplished simply by fitting 
the appropriate volume dependences to the results of the LQCD calculations for a range of 
d. From eq. (|9]), the binding momentum for any given d is 

K^"^ = Ko + ^ F^''\koL) + 0(e-2««V^) , (12) 

with the coefficients and kinematic factors in F^^\kqL) determined by the lattice calculation. 
Therefore, up to O (^e~'^'^°^ / Lj , the two free-parameters that remain to be determined are 

and which can be accomplished with a x^-minimization. In the case of having k^'^^ for 
only two d's, kq and can be solved for within the uncertainties of the LQCD calculations ^. 



V. CONCLUSIONS 

As recently stressed by Bour et al, the binding energy of a bound state depends upon its total 
momentum when subject to periodic boundary conditions in the spatial directions. Through 
the momentum modes excluded by the boundary conditions, these volume modifications of 
the binding energy depend upon the ratio of the spatial extent of the volume to the Lorentz- 
contracted size of the bound state, and also upon the masses of the constituents. We have 
extended the work of Bour et al from nonrelativistic quantum mechanics to quantum field 
theory and have pointed out that these features can be utilized in Lattice QCD calculations 
of hadronic bound states to (approximately exponentially) reduce the volume modifications 
of predicted binding energies. The standard Lattice QCD methodology that is used to 
determine the binding energy of a bound state is to measure the ground state energy of a 
system in a number of lattice volumes and then extrapolate to L = oo with a function of the 
form ~ e-'">^/L + (^e-^''°^/Ly Using combinations of boosted ground state energies, the 
volume dependence of the binding energy can be exponentially reduced in the nonrelativistic 
limit. For instance, the ground state energies of the lowest four boosted states in the lattice 
volume can be combined to reduce the volume modifications to the predicted binding energy 
to ~ e-2''"^/L + O (rfe-'^^^/L) where r/ < 1. 

In the specific case of the deuteron (neglecting its d-wave component), we have numer- 
ically explored what might be expected from future Lattice QCD calculations, and in par- 
ticular, examined the volume dependence of combinations of boosted ground state energies. 
We find that the deuteron binding energy can be extracted to high precision in reasonably 
modest volumes, reducing the volume modifications by more than an order of magnitude 



^ The volume extrapolation performed in Ref. [1] to determine the H-dibaryon binding energy, used the 
ground state energy obtained in two different lattice volumes and iteratively solved for and in 

^(0,0,0) 
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over those of the state at rest for L> 10 fm. It is also found that the volume modifications to 
the binding energy of the system with one unit of lattice momentum are significantly smaller 
than those of other low-momentum states, including the state at rest. It is clear that future 
Lattice QCD calculations that focus on extracting the properties and interactions of nuclei, 
including exotic systems such as the H-dibaryon, can greatly enhance the precision of their 
predicted binding energies by including systems with nonzero total momentum into their 
production. 

We would like to thank William Detmold and Dean Lee for useful discussions. ZD and 
MJS were supported in part by the DOE grant DE-FG03-97ER4014. 
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